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Abstract 

We investigate the in-medium masses of D and D mesons in the isospin-asymmetric nuclear 
matter at finite temperatures arising due to the interactions with the nucleons, the scalar isoscalar 
meson o", and the scalar iso-vector meson 5 within a SU(4) model. However, since the chiral 
symmetry is explicitly broken for the SU(4) case due to the large charm quark mass, we use 
the SU(4) symmetry here only to obtain the interactions of the D and D mesons with the light 
hadron sector, but use the observed values of the heavy hadron masses and empirical values of 
the decay constants. The in-medium masses of J /ip and the excited charmonium states ('0(3686) 
and i/;(3770)) are also calculated in the hot isospin asymmetric nuclear matter in the present 
investigation. These mass modifications arise due to the interaction of the charmonium states with 
the gluon condensates of QCD, simulated by a scalar dilaton field introduced to incorporate the 
broken scale invariance of QCD within the effective chiral model. The change in the mass of J /ip in 
the nuclear matter with the density is seen to be rather small, as has been shown in the literature 
by using various approaches, whereas, the masses of the excited states of charmonium (■i/'(3686) and 
'0(3770)) are seen to have considerable drop at high densities. The present study of the in-medium 
masses of D (D) mesons as well as of the charmonium states will be of relevance for the observables 
from the compressed baryonic matter, like the production and collective flow of the D (D) mesons, 
resulting from the asymmetric heavy ion collision experiments planned at the future facility of the 
FAIR, GSI. The mass modifications of D and D mesons as well as of the charmonium states in hot 
nuclear medium can modify the decay of the charmonium states (^ , Xc J /^) to DD pairs in the 
hot dense hadronic matter. The small attractive potentials observed for the D mesons may lead 
to formation of the D mesic nuclei. 
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INTRODUCTION 



The study of in-inediuin properties of hadrons is important for the understanding of strong 
interaction physics. The study of in-medium hadron properties has relevance in heavy-ion 
colhsion experiments as well as in nuclear astrophysics. There have been also extensive 
experimental efforts for the study of in-medium hadron properties by nuclear collision ex- 
periments. In these heavy-ion collision experiments, hot and dense matter is produced. By 
studying the experimental observables one can infer about how the hadron properties are 
modified in the medium. For example, the observed enhanced dilepton spectra l|-l3| could 
be a signature of medium modifications of the vector mesons [^-Isl . Similarly the properties 
of the kaons and antikaons have been studied experimentally by KaoS collaboration and the 
production of kaons and antikaons in the heavy-ion collisions and their collective flow are 
directly related to the medium modifications of their spectral functions p, |9Hl5[. The study 
of D and D mesons properties will be of direct relevance for the upcoming experiment at 
FAIR, GSI, where one expects to produce matter at high densities and moderate tempera- 
tures [l6| . At such high densities, the properties of the D and D mesons produced in these 
experiments are expected to be modified which should reflect in experimental observables 
like their production and propagation in the hot and dense medium. The reason for an ex- 
pected appreciable modifications of the D and D mesons is that D and D mesons contain a 
light quark (u,d) or light antiquark. This light quark or antiquark interacts with the nuclear 
medium and leads to the modifications of D and D properties. The experimental signature 
for this can be their production ratio and also in- medium J /ip suppression 17Hl9l| . In heavy- 
ion collision experiments of much higher collision energies, for example in RHIC or LHC, 
it is suggested that the J/-?/' suppression is because of the formation of quark-gluon plasma 



(QGP) |20|, I21J. However, in Ref. [22 



24j it is observed that the effect of hadron absorp- 



tion of J/ip is not negligible. In Ref. 25|, it was reported that the charmonium suppression 
observed in Pb + Pb collisions of NA50 experiment cannot be simply explained by nucleon 
absorption, but needs some additional density dependent suppression mechanism. It was 
suggested in these studies that the comover scattering |25l-l27] can explain the additional 
suppression of charmonium. An important difference between J /tp suppression pattern in 
comovers interaction model and in a deconfining scenario is that, in the former case, the 



to central collisions rather than in a 



anomalous suppression sets in smoothly from periphera 
sudden way when the deconfining threshold is reached 26|. The J /ip suppression in nuclear 
collisions at SPS energies has been studied in covariant transport approach HSD in Ref.(27J. 
The calculations show that the absorption of J/^'s by both nucleons and produced mesons 
can explain reasonably not only the total J /ip cross-section but also the transverse energy 
dependence of J /ip suppression measured in both proton-nucleus and nucleus collisions. In 
Ref. [28|, the cross section of J/ip dissociation by gluons is used to calculate the J/ip sup- 
pression in an equilibrating parton gas produced in high energy nuclear collisions. The large 
average momentum in the hot gluon gas enables gluons to break up the J /ip., while hadron 
matter at reasonable temperature does not provide sufficiently hard gluons. 

Due to the reduction in the masses of D and D mesons in the medium it is a possibility 
that excited charmonium states can decay to DD pairs [29i] instead of decaying to lowest 
charmonium state J /ip. Actually higher charmonium states are considered as major source 
of J /ip 3Cl|. Even at certain higher densities, it can become a possibility that the J /ip itself 
will decay to DD pairs. So this can be an explanation of the observed J /ip suppression by 
NA50 collaboration at 158 GeV/nucleon in the Pb-Pb collisions [20|. The excited states of 
charmonium also undergo mass drop in the nuclear medium 3l|]. The modifications of the 
in- medium masses of D mesons is large then the J /ip mass modification 32|, |33|. This is 
because the charmonium states are made up of a heavy charm quark and a charm antiquark. 
Within QCD sum rules, it is suggested that these heavy charmonium states interact with 
the nuclear medium through the gluon condensates. This is contrary to the interaction of 
the light vector mesons (p, w, 0), which interact with the nuclear medium through the quark 
condensates. This is because all the h eavy quark condensates can be related to the gluon 



condenstaes via heavy-quark expansion 3J]. Also in the nuclear medium there are no valence 



charm quark to leading order in density and any interaction with the medium is gluonic. 
The QCD sum rule approach 35| and leading order perturbative calculations 36[ to study 



the medium modifications of charmonium, show that the mass of J /ip is reduced slightly in 
the nuclear medium. In 3l|] , the mass modification of charmonium has been studied using 
leading order QCD formula and the linear density approximation for the gluon condensate in 
the nuclear medium. This shows a small drop for the J /ip mass at the nuclear matter density. 



but there is seen to be significant shift in the masses of the excited states of charmonium 
(^(3686) and ^(3770)). 

The in-medium modifications of D and D mesons have been studied using various ap- 
proaches. For example, in the QCD sum rule approach, it is suggested that the light quark or 
antiquark of D{D) mesons interacts with the light quark condensate leading to the medium 



modification of the D{D) meson masses 37|, |38|]. The quark meson coupling (QMC) model 



has also been used to study the D- meson properties [39|]. In the QMC model, the light 
quarks (u,d) and light antiquarks {u,d) confined in the nucleons and mesons interact via 
exchange of a scalar-isoscalar a meson as well a vector u meson. The nucleon has a large 
reduction of mass in the dense medium arising due to the interaction of the light quarks 
(u,d) with the a field. The drop in the mass of D mesons observed in the QMC model turns 
out to be similar to those calculated within the QCD sum rule approach. 

In the present investigation, we study the properties of the D and D mesons in the 
isospin-asymmetric hot nuclear matter. These modifications arise due to their interactions 
with the nucleons, the non-strange scalar isoscalar meson a and the scalar isovector meson 
6. We also study the medium modification of the masses of J/ip and excited charmonium 
states '?/'(3686) and '?/'(3770) in the nuclear medium due to the interaction with the gluon 
condensates using the leading order QCD formula. The gluon condensate in the nuclear 
medium is calculated from the medium modification of a scalar dilaton field introduced 



within a chiral SU(3) model 40| through a scale symmetry breaking term in the Lagrangian 
density leading to the QCD trace anomaly. In the chiral SU(3) model, the gluon condensate 
is related to the fourth power of the dilaton field x ^i^nd the changes in the dilaton field with 
the density are seen to be small. We study the isospin dependence of the in-medium masses 
of charmonium obtained from the dilaton field, x calculated for the asymmetric nuclear 
matter at finite temperatures. The medium modifications of the light hadrons (nucleons 
and scalar mesons) are described by using a chiral SU{3) model 40[. The model has been 
used to study finite nuclei, the nuclear matter properties, the in-medium properties of the 



vector mesons 



41 



in nuclear matter 



43 



42 1 as well as to investigate the optical potentials of kaons and antikaons 



44j | and in hyperonic matter in 45| . For the study of the properties 



D mesons in isospin-asymmetric medium at finite temperatures, the chiral SU(3) model is 



generalized to SU{A) flavor symmetry to obtain the interactions of D and D mesons with the 
light hadrons. Since the chiral symmetry is explicitly broken for the SU(4) case due to the 
large charm quark mass, we use the SU(4) symmetry here only to obtain the interactions of 
the D and D mesons with the light hadron sector, but use the observed values of the heavy 
hadron masses and empirical values of the decay constants. This has been in line with the 
philosophy followed in Ref. [46| where charmonium absorption in nuclear matter was studied 
using the SU(4) model to obtain the relevant interactions. However, the values of the heavy 
hadron masses and the coupling constants in Ref. J46| , were taken as the empirical values or 
as calculated from other theoretical models. The coupling constants were derived by using 
the relations from SU(4) symmetry, if neither the empirical values nor values calculated 
from other theoretical models were available [46|. The D meson prop erties in symmetric 
hot nuclear matter using SU(4) model have been studied in ref. 47|] and for the case of 



In a coupled channel approach 



asymmetric nuclear matter at zero temperature in 
for the study of D mesons, using a separable potential, it was shown that the resonance 
Ac(2593) is generated dynamically in the 1=0 channel J49[ analogous to A(1405) in the 
coupled channel approach for the KN interaction 50|. The approach has been generalized 
to study the spectral density of the D- mesons at finite temperatures and densities 5l|, 
taking into account the modifications of the nucleons in the medium. The results of this 
investigation seem to indicate a dominant increase in the width of the D-meson whereas 
there is only a very small change in the D-meson mass in the medium 511] • However, these 
calculations [49|,|5l|, assume the interaction to be SU(3) symmetric in u,d,c quarks and ignore 
channels with charmed hadrons with strangeness. A coupled channel approach for the study 
of D-mesons has been developed based on SU(4) symmetry [52] to construct the effective 
interaction between pseudoscalar mesons in a 16-plet with baryons in 20-plet representation 
through exchange of vector mesons and with KSFR condition [53]. This model 52|] has been 
modified in aspects like regularization method and has been used to study DN interactions 



in Ref. 



54| . This reproduces the resonance Ac(2593) in the 1=0 channel and in addition 



generates another resonance in the 1=1 channel at around 2770 MeV. These calculations have 
been generalized to finite temperatures 55(] accounting for the in-medium modifications of 
the nucleons in a Walecka type a — u model, to study the D and D properties 56| in the hot 



and dense hadronic matter. At the nuclear matter density and for zero temperature, these 
resonances (Ac(2593) and Sc(2770)) are generated 45 MeV and 40 MeV below their free 
space positions. However at finite temperature, e.g., at T = 100 MeV resonance positions 
shift to 2579 MeV and 2767 MeV for A^ (/ = 0) and E^ (/ = 1) respectively. Thus at finite 
temperature resonances are seen to move closer to their free space values. This is because 
of the reduction of pauli blocking factor arising due to the fact that fermi surface is smeared 
out with temperature. For D mesons in coupled channel approach a small repulsive mass 
shift is obtained. This will rule out of any possibility of charmed mesic nuclei 55| suggested 
in the QMC model 39|. But as we shall see in our investigation, we obtain a small attractive 



mass shift for D mesons which can give rise to the possibility of the formation of charmed 
mesic nuclei. The study of D meson self-energy in the nuclear matter is also helpful in 
understanding the properties of the charm and the hidden charm resonances in the nuclear 



matter 



571 ]. In coupled channel approach the charmed resonance Dso(2317) mainly couples 



to DK system, while the Do(2400) couples to D-k and DsK. The hidden charm resonance 
couples mostly to DD. Therefore any modification of D meson properties in the nuclear 
medium will affect the properties of these resonances. In Ref. 58|, |59|], the DN interactions 
at low energies have been studied using a meson exchange model {u and p) in close analogy 
with the meson exchange model for the KN interactions GOj, supplemented with a short- 
distance contribution from one-gluon exchange. The scattering lengths for the 1=0 and 1=1 
channels for the DN interactions arising from the one gluon exchange are seen to be very 
close to the values in Ref. 56|. Generalizing the SU(4) models with vector meson exchange 
potentials to SU(8) spin-flavor symmetry, that treats the heavy pseudoscalar and vector 
mesons on equal footing-as required by heavy-quark symmetry, in Ref. 6l|], the charmed 
baryon resonances that are generated dynamically have been studied within a unitary meson- 
baryon coupled-channel model. Some of the resonances in this model are identified with the 
recently observed baryon resonances. In the present investigation, the D{D) energies are 
modified due to a vectorial Weinberg- Tomozawa, scalar exchange terms {a, 6) as well as 



range terms 4J, |45|. The isospin asymmetric effects among D^ and D^ in the doublet, 



D= {D^, D^) as well as between D^ and D in the doublet, D = {D'^, D ) arise due to the 
scalar-isovector 6 meson, due to asymmetric contributions in the Weinberg- Tomozawa term. 



as well as in the range term 44 1. 



We organize the paper as follows. In section II, we give a brief introduction to the 
effective chiral SU{3) model used to study the isospin asymmetric nuclear matter at finite 
temperatures, and its extension to the SU{4) model to derive the interactions of the charmed 
mesons with the light hadrons. In section III, we present the dispersion relations for the D 
and D mesons to be solved to calculate their optical potentials in the hot and dense hadronic 
matter. In secton IV, we show how the in-medium masses of the charmonium states J/ip, 
'?/'(3686) and iplSTTO) in the present investigation, arise from the medium modification of the 
scalar dilaton field, introduced within the chiral model to incorporate broken scale invariance 
leading to QCD trace anomaly. Section V contains the results and discussions and finally, in 
section VI, we summarize the results of present investigation and discuss possible outlook. 

II. THE HADRONIC CHIRAL SU{3) x SU{3) MODEL 

We use a chiral SU{3) model for the study of the light hadrons in the present investigation 



|40| |. The model is based on nonlinear realization of chiral symmetry 62h64| and broken 



scale invariance 40l-l42l| . The effective hadronic chiral Lagrangian contains the following 
terms 

^ = ^kin + 2^ ^BW + ^vec + -^0 + ^SB (1) 

W=X,Y,V,Au 

In Eq.(l), Ckin is the kinetic energy term, Cbw is the baryon- meson interaction term in which 
the baryons-spin-0 meson interaction term generates the baryon masses. Cyec describes the 
dynamical mass generation of the vector mesons via couplings to the scalar mesons and 
contain additionally quartic self-interactions of the vector fields. Co contains the meson- 
meson interaction terms inducing the spontaneous breaking of chiral symmerty as well as a 
scale invariance breaking logarthimic potential. Csb describes the explicit chiral symmetry 
breaking. 

To study the hadron properties at finite temperature and densities in the present inves- 
tigation, we use the mean field approximation, where all the meson fields are treated as 
classical fields. In this approximation, only the scalar and the vector fields contribute to the 
baryon-meson interaction, Cbw since for all the other mesons, the expectation values are 



8 



zero. The interactions of the scalar mesons and vector mesons with the baryons are given 
as 



£ 



Bscal 



+ c 



Bvec 



Yl '^i \.'^*i + 5'a;i70W + 9iplOP + 5'0i7o0] A- 



(2) 



The interaction of the vector mesons, of the scalar fields and the interaction corresponding 
to the explicitly symmetry breaking in the mean field approximation are given as 



C: 



/ 2 2 2 2 2 / 2\ X 

- [m^u + nipp + m^0 j - 



xl 



+ ^4(c^' + 6u;V + p^ + 204), 



(3) 



C, = -\hx' (a^ + C^ + 5') + k, [a^ + e + 5' 



a 



4 ^4 



+ h[- + - + ^aH^ + C'\+hx(o'-5')C 



^4X^-^X^ln^ + ^X^ln( 



X^ d 
Xo 3 



{^'-s')C \rxy' 



and 



C 



SB 



X_ 
.Xo, 



TTT'lfnO- + V2mlfk T^mlf^ ( 



1 

V2' 



In dS]), rrii* is the effective mass of the baryon of species i, given as 



rrii 



-{9ai(y + 9ctC + 95t^) 



(4) 



(5) 



(6) 



The baryon-scalar meson interactions, as can be seen from equation (jH]), generate the baryon 
masses through the coupling of baryons to the non-strange a, strange C, scalar mesons and 
also to scalar-isovector meson 5. In analogy to the baryon-scalar meson coupling there 
exist two independent baryon-vector meson interaction terms corresponding to the F-type 
(antisymmetric) and D-type (symmetric) couplings. Here antisymmetric coupling is used 
because the universality principle J65| and vector meson dominance model suggest small 
symmetric coupling. Additionally, we choose the parameters 40|, |4^ so as to decouple the 
strange vector field 0^ ~ S7^s from the nucleon, corresponding to an ideal mixing between 
u and (j) mesons. A small deviation of the mixing angle from ideal mixing 66|-|68| has not 
been taken into account in the present investigation. 



9 



The concept of broken scale invariance leading to the trace anomaly in (massless) QCD, 
Qij^ = I'^^ C^^G^''", where GJ^^ is t^ gluon field strength tensor of QCD, is simulated in the 
effective Lagrangian at tree level 



C 



69| through the introduction of the scale breaking terms 

4\ A // T x/,.x3 



scalebreaHng ^X m^^^^J + ^X ^^\^Vdet(X)o 



x_ 

Xo 



(7) 



where J3 = det(X), with X as the multiplet for the scalar mesons. These scale breaking 
terms, in the mean field approximation, are given by the last two terms of the Lagrangian 
denstiy, Cq given by equation (jlj). The effect of these logarithmic terms is to break the scale 
invariance, which leads to the trace of the energy momentum tensor as 70 1 






(8) 



Hence the scalar gluon condensate of QCD (— (G"^,yC"^")) is simulated by a scalar dilaton 
field in the present hadronic model. 

The coupled equations of motion for the non-strange scalar field a, strange scalar field C,, 
scalar-isovector field 5 and dilaton field x, are derived from the Lagrangian density and are 
given as 



koxa - 4A;i [a' + C + 5')a- 2k2 [a' 

2 



3a6^) - 2hx(rC 



d 



2a 



3^ \a^-5^ 



+ (£\ mlU-J2g^,pt = 



(9) 



kox'C - 4A;i {a' + (^ + 6') ( - 4^^ - hx {cr' - 6' 



3 C \Xo, 



^mlfk - -J^^lf^ 



J29(^pt = Q 



(10) 



kox'5 - Ak, [a^ 
2 , J S \ 



+ 7;dx 



a' 



-52 



C' + 6^)6- 2k2 [6'' + 3a^6) + k^xK 



(11) 



koX {cr' + e + S') - k, {a' - 6') ( + x' 



1 + ln 






- TT^x'ln 



^oCo /Vxo/ 



2x 
Xo 



+ (4A;4 - d)x^ 
1 



mlU(r + f V2mlfk - -j^lh 1 C 



0(12) 



10 



In the above, pi* are the scalar densities for the baryons, given as 

^ _ r d?k m* f 1 1 \ 

^' ~ ^V (27r)3 E*{k) \e{E.*{k)-^..*)/T + i + ^(E-ik)+^.,*)/T + ij ^^'^) 

where, E,*{k) = [k"^ + mj*^)"'^/^, and, n* = /ij — (/^jo; — (/p^p — g^icf), are the single particle 
energy and the effective chemical potential for the baryon of species z, and, 7j=2 is the spin 



degeneracy factor [44 1. 



The above coupled equations of motion are solved to obtain the density and temperature 
dependent values of the scalar fields (a, C, and 5) and the dilaton field, Xi i^i the isospin 
asymmetric hot nuclear medium. As has already been mentioned, the value of the x is related 
to the scalar gluon condensate in the hot hadronic medium, and is used to compute the in- 
medium masses of charmonium states in the present investigation. The isospin asymmetry 
in the medium is introduced through the scalar-isovector field 5 and therefore the dilaton 
field obtained after solving the above equations is also dependent on the isospin asymmetry 
parameter, rj defined as r/ = (p„ — pp)/{2pB), where p„ and Pp are the number densities 
of the neutron and the proton and pb is the baryon density. In the present investigation, 
we study the effect of isospin asymmetry of the medium on the masses of the charmonium 
states J/^,?/'(3686) and ^(3770). 

The comparison of the trace of the energy momentum tensor arising from the trace 
anomaly of QCD with that of the present chiral model gives the relation of the dilaton field 



to the scalar gluon condensate. We have, in the limit of massless quarks 71 1 



^, = ^p^G;^G^n ^ -(1 - d)x' (14) 

The parameter d originates from the second logarithmic term of equation (jTj). To get an 
insight into the value of the parameter d, we recall that the QCD /3 function at one loop 
level, for Nc colors and Nf flavors is given by 

^-<^)-^ 0-1^) -«<«=) (-) 

In the above equation, the first term in the parentheses arises from the (antiscreening) self- 
interaction of the gluons and the second term, proportional to Nf, arises from the (screening) 
contribution of quark pairs. Equations flT^ and flTSl) suggest the value of d to be 6/33 for 
three flavors and three colors, and for the case of three colors and two flavors, the value 

11 



of d turns out to be 4/33, to be consistent with the one loop estimate of QCD /3 function. 



These values give the order of magnitude about which the parameter d can be taken 70| . 
since one cannot rely on the one-loop estimate for /3qcd(5')- In the present investigation of 
the in-medium properties of the charmonium states due to the medium modification of the 



dilaton field within chiral SU{?>) model, we use the value of (i=0.064 [48|]. This parameter, 
along with the other parameters corresponding to the scalar Lagrangian density, Cq given 
by dl]), are fitted so as to ensure extrema in the vacuum for the a, C, and x field equations, 
to reproduce the vacuum masses of the rj and rj' mesons, the mass of the a meson around 



500 MeV, and pressure, p(po)=0, with po as the nuclear matter saturation density J40|, l48 |. 
The trace of the energy- momentum tensor in QCD, using the one loop beta function 
given by equation flTSl) . for A^c=3 and Nf=3, is given as, 

0, = --fG^G^'^^ (16) 



Using equations (1141) and (1161) . we can write 



-IG^^G^^'^) = -{I - d)x' (17) 



~ -i^-- I 9 

We thus see from the equation (ITTj) that the scalar gluon condensate i^G'^ G^'^"'\ is propor- 
tional to the fourth power of the dilaton field, X) in the chiral SU(3) model. As mentioned 
earlier, the in-medium masses of charmonium states are modified due to the gluon conden- 
sates. Therefore, we need to know the change in the gluon condensate with density and 
temperature of the asymmetric nuclear medium, which is calculated from the modification 
of the X field, by using equation (TT7I) . 

III. D AND D MESONS IN HOT ASYMMETRIC NUCLEAR MATTER 

In this section we study the D and D mesons properties in isospin-asymmetric nuclear 
matter at finite temperatures. The medium modifications of the D and D mesons arise due 
to their interactions with the nucleons and the scalar mesons and the interaction Lagrangian 



density is given as J48| 



Cvu - -g^ 



12 



pYp - nYn) (/^°(9^5°) - (9^/^0)5°) - (D+id^D-) - {d,D+)D- 



+ 



m 



D 



Vd 
1 



/ 



D 



[a + V2Cc)P°I^° + {D-D+)) + (5(5°D°) - {D-D+)) 



a + v^Cc)((5^5°)(9^D°) + {d^D~){d^D- 



+ 5((a^L)°)(a'^D°) - (9^D-)(9^D^ 



di 



+ 



2/|> 
4^ 






+ (pp - nn){{d^D^){d^D'^)) - {d^D-){d^D+)) 



(18) 



In Eq.l TTSl) . the first term is the vectorial Weinberg Tomozawa interaction term, obtained 
from the kinetic term of Eq.(l). The second term is obtained from the exphcit symmetry 
breaking term and leads to the attractive interactions for both the D and D mesons in the 
medium. The next three terms of above Lagrangian density (~ {d^D){d^D)) are known 
as the range terms. The first range term (with coefficient ( — j-)) is obtained from the 
kinetic energy term of the pseudoscalar mesons. The second and third range terms di and 
d2 are written for the DN interactions in analogy with those written for KN interactions 
in 45|. It might be noted here that the interaction of the pseudoscalar mesons with the 
vector mesons, in addition to the pseudoscalar meson-nucleon vectorial interaction, leads to 
a double counting in the linear realization of chiral effective theories. Further, in the non- 
linear realization, such an interaction does not arise in the leading or subleading order, but 
only as a higher order contribution 72|. Hence the vector meson-pseudoscalar interactions 
will not be taken into account in the present investigation. 

The dispersion relations for the D and D mesons are obtained by the Fourier transfor- 
mations of equations of motion. These are given as 



- w^ + A^2 + m^ - n (w, \k\) = 



(19) 



where, mo is the vacuum mass of the D{D) meson taken as 1869 MeV and 1864.5 MeV for 
the D and D mesons respectively. 11 (to, \k\] denotes the self-energy of the D (D) mesons 
in the medium. 

The self-energy 11 (to, \k\] for the D meson doublet {D'^, D^) arising from the interaction 



13 



of Eq. fjTSjl is given as 

Uioj,\k\) 



4/1. 



m 



S{pp + Pn) ± {Pp - Pn) 



u 



+ ^{cr' + V2Cc'±S') 



2/ 



D 



+ 



+ 



Id ^Jd 



d2 



p% + p\)±{p%-p\)) iu'-k'), 



(20) 



where the ± signs refer to the D^ and D~^ mesons, respectively, and cr'(cr — o"o), 
C (Cc — Cco), and 6' {= 6 — 6q) are the fluctuations of the scalar isoscalar fields a and ( and 
the scalar-isoscalar field 6 from their vacuum expectation values. The vacuum expectation 
value of 6 is zero {6o = 0), since a nonzero value for it will break the isospin-symmetry of the 
vacuum. (We neglect here the small isospin breaking effect arising from the mass and charge 
difference of the up and down quarks.) We might note here that the interaction of the scalar 
quark condensate Cc (being made upof heavy charmed quarks and antiquarks) leads to very 
small modifications of the masses J73| . So we will not consider the medium fluctuations of 
(c- In the present investigation, we take the value of the D meson decay constant, fr, as 135 
MeV [7^. Within the present model, the medium modification to the D{D) mesons due to 
the scalar interaction depends only on the fluctuations of the scalar a and 6 fields in the 
asymmetric hot nuclear medium which are determined by solving the coupled equations (j9]), 
( ITOj) . ( ITTj) and ( 1T21) . As the scalar term of ( 120|) does not contain any free parameters, with 
the assumption that the fluctuation of the charm condensate has a negligible effect on the 
masses of the D{D) mesons, there are no uncertainties in the mass shift due to the scalar 
interaction in the present investigation, once the value for f^ is chosen. In Eq. fl20|) . pp and 
pn are the number densities of protons and neutrons given by 

d^k f 1 1 \ 



for z=p and n, and p^ and p^ are their scalar densities, as given by equation flT^ . 
Similarly, for the D meson doublet (D^,D~j, the self-energy is calculated as 

1 
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(22) 



where the ± signs refer to the D^ and D mesons, respectively. The optical potentials of 
the D and D mesons are obtained using the expression 



U{u,k) = uJ{k)-^Jk^ + m 



(23) 



where mo is the vacuum mass for the D{D) meson and uj{k) is the momentum-dependent 
energy of the D{D) meson. 

IV. CHARMONIUM MASSES IN HOT ASYMMETRIC NUCLEAR MATTER 



In this section, we investigate the masses of charmonium states J/ip, '?/'(3686) and '?/'(3770), 
in isospin asymmetric hot nuclear matter. From the QCD sum rule calculations, the mass 

nn 

shift of the charmonium states in the medium is due to the gluon condensates [31l, |37[ . For 
heavy quark systems, there are two independent lowest dimension operators: the scalar gluon 
condensate U^G'^^^G^^^"-^) and the condensate of the twist 2 gluon operator {{-^G'^^^G^""')) . 
These operators can be rewritten in terms of the color electric and color magnetic fields, 
{—E"^) and (—5^). Additionally, since the Wilson coefficients for the operator {—B"^) 
vanish in the non-relativistic limit, the only contribution from the gluon condensates is 
proportional to {—E"^), similar to the second order Stark effect. Hence, the mass shift of 
the charmonium states arises due to the change in the operator {—E"^) in the medium from 
its vacuum value 3l|. In the leading order mass shift formula derived in the large charm 



mass limit 



the shift in the mass of the charmonium state is given as 3l|] 

k 



Am4e) = -- J dk \^^\ 



k'^/rric + e V \ TT 



—E^ 

TT 



(24) 



In the above, rric is the mass of the charm quark, taken as 1.95 GeV 3l|, ''^v ^^ ^^^ ^^'^~ 
uum mass of the charmonium state and e = 2r7ic — rn^. ip{k) is the wave function of the 
charmonium state in the momentum space, normalized as J j:;^\4'{k)\'^ = 1 75|. At finite 



15 



densities, in the linear density approximation, the change in the value of {—E"^), from its 
vacuum value, is given as 



77 / \ TT 
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TT 



PB 

n2Mn 



(25) 



and the mass shift in the charmonium states reduces to 31 1 
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In the above, i^E ) is the expectation value of i^E ) with respect to the nucleon. 

The expectation value of the scalar gluon condensate can be expressed in terms of the 
color electric field and the color magnetic field as 
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(27) 



In the non-relativistic limit, as already mentioned, the contribution from the magnetic field 
vanishes and hence, we can write, 
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« 7712 
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(28) 



Using equations ( TT7|) . 1^^ and ( 12S]) . we obtain the expression for the mass shift in the 
charmonium in the hot and dense nuclear medium, which arises from the change in the 
dilaton field in the present investigation, as 



4 r 

Am^e) = —{1 - d) J dk' 



M{k) 



/ 4 4 

[X -Xo 



(29) 



dk k^/mc + e 

In the above, x ^^^ Xo ^i-re the values of the dilaton field in the nuclear medium and in the 
vacuum respectively. 



In the present investigation, t 



le wave functions for the charmonium states are taken to 



be Gaussian and are given as 33 1 



la2^2^l+l 



^N,i = Normalization x Yi"'{9,^){l3\-^)^^exp-^'^^''^ Lj^^^ (/^V^) 



(30) 



where /3^ = Mu/h characterizes the strength of the harmonic potential, M = mc/2 is the 
reduced mass of the charm quark and charm anti-quark system, and Vl{z) is the associ- 
ated Laguerre Polynomial. As in Ref. 3l|, the oscillator constant /3 is determined from 
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the mean squared radii (r^) as 0.46^ fm^, 0.96^ fm^ and 1 fm^ for the charmonium states 
J/'?/'(3097), ^(3686) and ip{3770), respectively. This gives the value for the parameter (3 as 
0.51 GeV, 0.38 GeV and 0.37 GeV for J/V^(3097), ^(3686 and V^(3770), assuming that these 
charmonium states are in the IS, 2S and ID states respectively. Knowing the wave functions 
of the charmonium states and calculating the medium modification of the dilaton field in 
the hot nuclear matter, we obtain the mass shift of the charmonium states, J/ip, '?/'(3686) 
and '?/'(3770) respectively. In the next section we shall present the results of the present 
investigation of these in-medium charmonium masses in hot asymmetric nuclear matter. 

V. RESULTS AND DISCUSSIONS 

In this section, we present the results of our investigation for the in-medium masses of 
D and D mesons as well as of the charmonium states J/'?/' (3097), '?/'(3686) and '?/'(3770), in 
isospin asymmetric nuclear matter at finite temperatures. We have generalized the chiral 
SU{3) model to SU{4) to include the interactions of the charmed mesons. The values 
of the parameters used in the present investigation, are : fco = 2.54, /ci = 1.35, A;2 =- 
4.78, ks = —2.77, k^ = —0.22 and d = 0.064, which are the parameters occurring in the scalar 
meson interactions defined in equation (jl]). The vacuum values of the scalar isoscalar fields, 
a and ( and the dilaton field x ^^e —93.3 MeV, —106.6 MeV and 409.8 MeV respectively. 
The values, g^j^ = 10.6 and gi^^ = —0.47 are determined by fitting to vacuum baryon masses. 
The other parameters fitted to the asymmetric nuclear matter saturation properties in the 
mean-field approximation are: g^^ = 13.3, gpp = 5.5, g^ = 79.7, gsp = 2.5, m^ = 1024.5 
MeV, mcr = 466.5 MeV and ms = 899.5 MeV. The nuclear matter saturation density used 
in the present investigation is 0.15 fm^^. The coefficients di and d2, calculated from the 
empirical values of the KN scattering lengths for J = and 1 = 1 channels, are 2.56/mK 
and 0.73/mi^, respectively J45|. 

In isospin asymmetric nuclear medium, the properties of the D mesons {D^, D^) and D 
mesons {D^,D^), due to their interactions with the hot hadronic medium, undergo medium 
modifications. These modifications arise due to the interactions with the nucleons (through 
the Weinberg- Tomozawa vectorial interaction as well as through the range terms) and the 
scalar exchange terms. The modifications of the scalar mean fields modify the masses of the 
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nucleons in the hot and dense hadronic medium. Before going into the details of how the D 
and D mesons properties are modified at finite temperatures in the dense nuclear medium, 
let us see how the scalar fields are modified at finite temperatures in the nuclear medium. 
In figures [H E], [3] and HI we show the variation of scalar fields a, C, and scalar-isovector 
field 5 and the dilaton field x, with temperature, for both zero and finite densities, and for 
selected values of the isospin asymmetry parameter, r^ = 0, 0.1, 0.3 and 0.5. At zero baryon 
density, we observe that the magnitudes of the scalar fields a and C, decrease with increase 
in temperature. However, the drop in their magnitudes with temperature is negligible upto 
a temperature of about 100 MeV (that is, they remain very close to their vacuum values). 
The changes in the magnitudes of the a and C, fields are 2.5 MeV and 0.8 MeV respectively, 
when the temperature changes from 100 MeV to 150 MeV, above which the drop increases. 
These values change to about 10 MeV and 3 MeV for the a and C, fields respectively, for 
a change in temperature of 100 MeV to 175 MeV. At zero baryon density, it is observed 
that the value of the dilaton field remains almost constant upto a temperature of about 130 
MeV above which it is seen to drop with increase in temperature. However, the drop in the 
dilaton field is seen to be very small. The value of the dilaton field is seen to change from 
409.8 MeV at T = to about 409.7 MeV and 409.3 MeV at T = 150 MeV and T = 175 
MeV, respectively. The thermal distribution functions have an effect of increasing the scalar 
densities at zero baryon density, i.e. /i* = 0, as can be seen from the expression of the scalar 
densities, given by equation (IT3|l . This effect seems to be negligible upto a temperature 
of about 130 MeV. This leads to a decrease in the magnitudes of the scalar fields, a and 
C^. This behaviour of the scalar fields is reflected in the value of Xi which is solved from 
the coupled equations of motion of the scalar fields, given by equations ([9]), flTO|) . flTTj) and 
(TT2|) . as a drop as we increase the temperature above a temperature of about 130 MeV. The 
scalar densities attaining nonzero values at high temperatures, even at zero baryon density, 
indicates the presence of baryon-antibaryon pairs in the thermal bath and has already been 



observed in literature 



m 



TJ\. This leads to the baryon masses to be different from their 
vacuum masses above this temperature, arising from modifications of the scalar fields a and 

C- 

For finite density situation, the behaviour of the scalar fields with temperature is seen to 
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be very different from the zero density case, as can be seen from the subplots (b), (c) and 
(d) of figures [1], [2] and H] for a, ( and x fields respectively, where the fields are plotted as 
functions of temperature for densities po, 2po and 4po respectively. At finite densities, one 
observes first a rise and then a decrease of the scalar fields a, ( and x with temperature. 
For example, at pb = Po, and for the value of the isospin asymmetry parameter, 77 = the 
scalar fields a, ( and x increase with temperature upto a temperature of about 145 MeV and 
above this temperature, they both start decreasing. At rj = 0.5 the value of temperature 
upto which these scalar fields increase becomes about 120 MeV. For ps = 4po and rj = the 
scalar fields a, ( and x increase upto a temperature of about 160 MeV. At rj = 0.5 this value 
of temperature is lowered to about 120 MeV for a, ( and x fields respectively. This observed 
rise in the magnitudes of a and ( fields with temperature leads to an increase in the mass of 
nucleons with temperature at finite densities. The reason for the different behaviour of the 
scalar fields (a and () at zero and finite densities can be understood in the following manner 



42]. As has already been mentioned, the thermal distribution functions in flT3ll have an 
effect of increasing the scalar densities at zero baryon density, i.e., for p* = 0. However, at 
finite densities, i.e. for nonzero values of the effective chemical potential, p*, for increasing 
temperature, there are contributions also from higher momenta, thereby, increasing the 
denominator of the integrand on the right hand side of equation flT3|) . This leads to a 
decrease in the scalar density. The competing effects of the thermal distribution functions 
and the contributions from the higher moments states give rise to observed behaviour of the 
scalar density and hence of the a and ( fields with temperature at finite baryon densities 



42 1 . This kind of behaviour of the scalar a field on temperature at finite densities has also 

n 

been observed in the Walecka model by Li and Ko [8]. The behaviour of the scalar fields a 
and ( with the temperature is reflected in the behaviour of x field, since it is solved from 
the coupled equations of the scalar fields. In figure [HI we observe that the value of the scalar 
isovector field, 6 is zero at zero baryon density, since there is no isospin asymmetry at zero 
density. At finite baryon densities the magnitude of the scalar isovector field, S decreases 
with increase in the temperature of the nuclear medium. However, for given temperature 
as we move to higher densities the magnitude of the S field is seen to increase which means 
that the medium is more asymmetric at higher densities. In the isospin asymmetric nuclear 
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medium, the scalar isovector 5 meson attains a nonzero expectation value. However, the 
magnitudes of the 5 field at a given baryon density and temperature of the medium, which 
are solved from the coupled equations, ([9]) to ( IT2|) . for the scalar fields, turn out to be small 
(< 4 — 5Mey),as can be seen from figure [3l This causes changes in the a and Q fields, as 
can be seen from figures [1] and [21 also to be small in the isospin asymmetric nuclear matter, 
since these solved from the equations (19]) and flTO|) . along with the equations for 5 and Xi 
given by (H]) and ([H]). 

In figures |5] and El we show the variation of the energy of D and D mesons due to 
Weinberg- Tomozawa term, the scalar exchange term arising from the explicit symmtery 
breaking term and the range terms, as functions of temperature and for the densities 
Pb = 0,Pq and 4po- These are plotted asymmetric nuclear matter with the value of the 
isospin asymmetry parameter as 0.5, and in the same figure, the results arising from these 
contributions for the isospin symmetric case (^=0) are shown as the dotted lines. At zero 
density, the contributions from the Weinberg- Tomozawa term is zero, since the densities of 
the proton and neutron, pp and p„ are zero at zero density. The scalar fields attaining values 
different from their vacuum values at zero density at values of temperature higher than 
about 100 MeV, as already observed in figures [1] and |2l leads to a decrease in the D^ and 
D^ mesons due to the scalar meson exchange term and an increase in the contribution due 
to the range term, as can be seen in the subplots (a) and (b) of figure El These two terms 
seem to almost cancel with each other, resulting in negligible modifications of D^ and D° 
mesons. The drop in the masses of D~^ and D^ are seen to be about 0.75 MeV and 0.75 MeV 
respectively at temperature T=150 MeV and about 5.4 MeV and 5.35 MeV at a temperature 
of 175 MeV. At finite densities, the contribution to the energies of the D^ and D^ mesons 
arising from the Weinberg- Tomozawa term as a function of temperature, is constant, since 
this term depends on the densities of proton and neutron only, which for a given baryon 
density and a given value of the isospin asymmetry parameter, are constant. At ps = po and 
T=0, the Weinberg- Tomozawa term gives a drop of about 23 MeV in the masses of D^ and 
D^ mesons from their vacuum values in isospin symmetric matter and about 31 MeV and 16 
MeV for D~^ and D° mesons in the isospin asymmetric matter with ?7=0.5. The mass drop 
of D+ meson of 23 MeV in symmetric nuclear matter at zero temperature may be compared 
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with the value of about 43 MeV, which is the drop in D^ mass in the Born approximation 
of the calculation of the coupled channel approach of Ref. [54|, when only the Weinberg 
Tomozawa interaction is included. In the presence of the scalar exchange term as well, the 
drop in the D^ mass in the present investigation is seen to be about 144 MeV, which may 
be compared to the drop in the mass of the D^ meson with Weinberg- Tomozawa as well as 
scalar interaction of about 60 MeV in Ref. [54]. In the present investigation, we observe 



the scalar interaction to be much more dominant than in Ref. [5J]. The contributions to 
the energy of D^ and D^ mesons due to the explicit symmetry breaking term (scalar me- 
son exchange term) are observed to increase with temperature upto a particular value of 
temperature, above which it is seen to decrease with temperature. This is a reflection of 
the behaviour of the scalar fields a and C, shown in figures [1] and [2], whose magnitudes are 
seen to increase with temperature upto a particular value, above which there is seen to be a 
decrease in their magnitudes. The range term is seen to give an increase in the D^ and D^ 
masses and at pB = po, there is seen to be a rise of about 71 MeV at zero temperature for 
?7=0. This value is modified to about 67 MeV at T=150 MeV and about 70 MeV at T=175 
MeV. For a density of ps = 4po and at zero temperature, the drop in the D-meson masses 
in the symmetric nuclear matter is seen to be about 27 MeV due to the range term. For 
this density, the contribution of the di and d2 terms dominate over the first range term thus 
leading to an overall drop of the D-meson masses in the medium of about 347 MeV in the 
symmetric nuclear matter at zero temperature. 

It may be noted that the difference in the in-medium energy of D^ and D~^ mesons arises 
from Weinberg- Tomozawa term only (see equations ( 120|) and ( 122|) ) whereas the the scalar 
meson exchange term and the range terms have identical contributions to the energy of D~^ 
and D~ mesons. The repulsive contribution for the D~ and D^ mesons for the Weinberg- 
Tomozawa term as compared to the attractive contribution for the D^ and D^ mesons, gives 
a rise in the masses of the D mesons (plotted in figure E]) as compared to the masses of the D 
mesons shown in figure [5l At zero density, the contribution from the Weinberg- Tomozawa 
term is zero and hence the modifications for the D~^ and D~ are identical, both having 
negligible changes in their masses due to the cancelling effects of the attractive scalar term 
and the repulsive range terms. At nuclear matter saturation density, when one takes into 
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account only the contribution from the Weinberg- Tomozawa term, there is seen to be an 
increase in the mass of the D~ meson of about 24 MeV. However, when we consider the 
contributions due to all the individual terms, the attractive scalar interaction is seen to 
dominate over the repulsive Weinberg- Tomozawa and the range terms, therby giving a drop 
of the D~ mass of about 27 MeV in the symmetric nuclear matter for pb = Po- However, 
for higher densities, the di and d2 terms of the range term, which are both attractive, give 
a further drop of the mass of the D~ meson and for density, pb = 4po, the overall drop is 
seen to be about 162 MeV for ri=0. 

We study the density dependence of D and D masses at finite temperatures at selected 
values of the isospin asymmetric paramater, t] and compare the results with the zero tem- 
perature case 481]. The isospin symmetric part of the Weinberg- Tomazawa term gives a 
drop of the D mass, as can be seen from the expressions of the self-energies of the D and 
D mesons given by equations ( !20l) and ( !22l) . However, the isospin asymmetric part of this 
term is seen to give a mass splitting for the D^ and D^, given by the second term of the 
Weinberg- Tomazawa term, giving a further drop of the D^ mass, whereas the asymmetry 
reduces the drop of the mass of D^. For the D mesons, the isospin symmetric part of the 
Weinberg Tomazawa term gives an increase in the mass and the isospin asymmetric con- 
tribution of this term gives a further rise in the mass of the D~ mass, whereas it reduces 
the increase of the D'^ mass in the asymmetric nuclear medium. For both D and D mesons 
in the symmetric nuclear matter, the scalar meson interaction is attractive and identical 
(except for a small difference due to the difference in the vacuum masses of D"^ and D^{D^) 
mesons). One might observe from the expressions of the D{D) self energies, given by equa- 
tions (|20|) and fl22|) that the non-zero value of 6 meson arising due to the isospin asymmtery 
in the medium gives a drop in the masses of D^{D^) in the asymmteric nuclear matter, 
whereas this interaction is repulsive for D^{D^) mesons. The contributions to the D and D 
masses due to the range terms are given by the last three terms of the self energies given by 
equations ( l20l) and ( l22l) . The first of these range terms is repulsive, whereas the second and 
third terms are attractive, when isospin asymmtery is not taken into account. However, due 
to a nonzero value of the 5 field arising from isospin asymmtery in the medium, the 5 term 
of the first range term leads an increase in the masses of the D^{D^) and a drop in the 
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FIG. 1: (Color online) The scalar-isoscalar field a plotted as a function of temperature at a given 
baryon density {pB = 0,po,2pQ and 4po), for different values of the isospin asymmetry parameter, 
V- 

masses of D^{D^) mesons. The second of the range terms (the di term) is attractive and 
gives identical mass drops for D^ and D^ in the D doublet as well as for D^ and D^ in the D 
doublet. This term is proportional to {p^^+p^), which turns out to be different for the isospin 
asymmetric case as compared to the isospin symmetric nuclear matter, due to the presence 
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FIG. 2: (Color online) The scalar-isoscalar field ^ plotted as a function of temperature at a given 
baryon density {pB = 0,po,2pQ and 4po), for different values of the isospin asymmetry parameter, 
77. 

of the S meson. This is because the equations of motion for the scalar fields for the two 
situations (with/without 6 mesons) give different values for the mean field, a (~ (ps^ + Ps^))- 
The last term of the range term (the ^2 term) has a negative contribution for the energies 
of D^ and D^ mesons as well as for D^ and D*^ mesons for the isospin symmetric matter. 
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FIG. 3: (Color online) The scalar-isovector field b plotted as a function of temperature at a given 
baryon density (p^ = 0,po,2pQ and 4po), for different values of the isospin asymmetry parameter, 
77. 

The isospin asymmetric part arising from the {{ps"' — Ps^)) term of the d2 term has a further 
drop in the masses for D^ mesons, whereas it increases the masses of the D^ and D'^ mesons 
from their isospin symmetric values. In Fig. [TJ we show the variation of the energy of the D 
mesons {D^,D^) at zero momentum with baryon density pb for different values of isospin 
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FIG. 4: (Color online) The dilaton field x plotted as a function of temperature at a given baryon 
density {ps = 0,po,2pQ and 4po), for different values of the isospin asymmetry parameter, 77. 

asymmetry parameter rj and with the values of temperature as T = 0, 100, 150 MeV. The 
isospin-asymmetry in the medium is seen to give an increase in the D^ mass and a drop in 
the D^ mass as compared to the isospin symmetric (77 = 0) case. This is observed both for 
zero temperature 48| and finite temperature cases. At nuclear matter saturation density, 
Pb = P01 the drop in the mass of D^ meson from its vacuum value (1869 MeV) is 78 MeV 
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FIG. 5: (Color online) The energies of D^ meson ((a),(c) and (e)) and of D^ meson ((b), (d) and 
(f)), at momentum k = 0, versus temperature, T, for different values of the isospin asymmetry 
parameter {rj = and 0.5) and for given values of density (ps = 0, Po and 4/9o)- The values of 
parameters di and d2 are calculated from KN scattering lengths in I = and 1 = 1 channels. 
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FIG. 6: (Color online) The energies of D^ meson ((a),(c) and (e)) and of D^ meson ((b), (d) and 
(f)), at momentum A: = 0, versus temperature, T, , for different values of the isospin asymmetry 
parameter {rj = and 0.5) and for given values of density (ps = 0, Po and 4po)- The values of 
parameters di and d2 are calculated from KN scattering lengths in I = and 1 = 1 channels. 
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for zero temperature case in isospin symmetric medium. At a density of 4po, this drop in 
the mass of D^ meson is seen to be about 347 MeV for ?7=0. At finite temperatures, the 
drop in the mass of D^ meson for a given value of isospin asymmetry decreases as compared 
to the zero temperature case. For example, at nuclear saturation density, po? the drop in 
the mass of D^ meson turns out to be 72, 65 and 62 MeV at a temperature of T = 50, 
100 and 150 MeV respectively for the isospin symmetric matter. At a higher density of 
Pb = 4:Pq, zero temperature value for the D~^ mass drop of about 347 MeV is modified to 
336, 313 and 294 MeV for T = 50, 100 and 150 MeV respectively. Thus the masses of D 
mesons at finite temperatures and finite densities are observed to be larger than the values 
at zero temperature case. This is because of the increase in the magnitudes of the scalar 
fields a and ( with temperature at finite densities as mentioned earlier. The same behaviour 
remains for the isospin asymmetric matter. This behaviour of the nucleons and hence of the 
D-mesons with temperature was also observed earlier for symmetric nuclear matter at finite 
temperatures within the chiral effective model J47J. The drop in the mass of -D^ meson is 
seen to be larger as we increase the value of the isospin asymmetry parameter. We observe 
that as we change rj from to 0.5, then the drop in the mass of D~^ meson is 95 MeV and 
384 MeV at densities of po and 4po respectively, for the zero temperature case. At T = 50 
MeV these values change to 89 MeV at po and 377 MeV at a density of 4po. For a baryon 
density, ps = Po(4po), the drop in the D+ mass is 83 (363) MeV at T = 100 MeV, and 84 
(359) MeV at T=150 MeV. We thus observe that for a baryon density of po, for ?7=0.5, the 
drop in the D~^ mass from its value for ri=0 is 17, 18 and 22 MeV for T=0, 100 and 150 MeV 
respectively, and for ps = 4po, these values are modifed to 37, 50 and 65 MeV respectively. 
Thus we observe that for a given value of density, as we move from rj = to rj = 0.5 the 
drop in mass of D'^ mesons is larger at higher temperatures. 

The mass of the D^ meson drops with density as can be seen from figure [7l The drop in 
the mass of D^ meson at density po from its vacuum value (1864.5 MeV) is 78, 65, 62 MeV 
for temperature T = 0, 100, 150 MeV respectively at r^ = 0. At density 4po and isospin- 
asymmetry parameter rj = 0, these values become 347, 313 and 294 MeV for temperatures 
T=0,100 and 150 MeV respectively. For the D^ meson, there is seen to be an increase in the 
mass as we move from isospin symmetric to isospin asymmetric medium. For example, at 
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FIG. 7: (Color online) The energies of D^ meson ((a),(c) and (e)) and of D^ meson ((b), (d) and 
(f)), at momentum k = 0, versus the baryon density (in units of nuclear saturation density), Pb/po, 
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FIG. 8: (Color online) The energies of D^ meson ((a),(c) and (e)) and of D^ meson ((b), (d) and 
(f)), at momentum k = 0, versus the baryon density, expressed in units of nuclear matter saturation 
density, pb/po, for different values of the isospin asymmetry parameter (rj = 0, 0.1, 0.3, 0.5) and for 
a given temperature (T = 0, 100 MeV and 150 MeV). The parameters di and d2 are determined 
from KN scattering lengths in I = and 1=1 channels. 
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zero temperature and baryon density equal to po and 4po, the rise in the masses of the Z)° 
meson are 21 and 91 MeV respectively as we move from isospin-symmetric medium (77 = 0) 
to isospin-asymmetric medium [r] = 0.5). At a temperature of 100 MeV, these values become 
17 MeV at po and 72 MeV at 4po. For T=150 MeV, these values become 9 MeV at po and 
44 MeV at 4po. Thus for D^ mesons, the rise in the mass, is seen to be lowered at higher 
temperatures as we move from rj = to rj = 0.5. The strong isospin dependence of the D~^ 
and D^ meson masses should show up in observables such as their production as well as flow 
in asymmetric heavy-ion collisions planned at the future facility at FAIR, GSI. 

FigJH] shows the results for the density dependence of the energies of the D mesons at 
zero momentum at values of the temperature, T = 0, 100, 150 MeV. There is seen to be 
a drop of the masses of both the D~ and D^ with density. This is due to the dominance 
of the attractive scalar exchange contribution as well as the range terms (which becomes 
attractive above a density of about 2-2.5 times the nuclear matter saturation density) over 
the repulsive Weinberg- Tomozawa interaction [48] • It is observed that the drop in the mass 
of D^ and D^ in isospin symmetric nuclear matter is 27.2 MeV at po and 162 MeV at 4po, 
at zero temperature, from their vacuum values. As we go to higher temperatures, the drop 
in the masses of D~ and -D° mesons decreases. For example, at r^ = and pb = Po the 
drop in the mass of D~ meson is 20.9, 14.3 and 10.8 MeV at a temperature of 50, 100 
and 150 MeV respectively. For temperature, T=0,50,100 and 150 MeV, at ps = Po, the 
values of the drop in the D'^ mass are seen to be 27.2, 21.4, 14.6 and 15.5 MeV respectively, 
and the drop in the D^ mass are modifed from about 27, 21, 14.2 and 11 MeV to 23, 19, 
14 and 19 MeV, when rj is changed from to 0.5. The masses of D~ and -D° mesons are 
observed to have negligible dependence on the isospin asymmetry upto a density of pb = Po- 
However, at high densities there is seen to be appreciable dependence of these masses on the 
parameter, rj. As we change 77 from to 0.5, at pb = 4po, the drop in the mass of the D~ 
meson is modified from 162 MeV to 139 MeV at zero temperature. At higher temperatures, 
T= 50, 100 and 150 MeV and at the density of 4po, when we change rj from to 0.5, the 
values of the D~ mass drop are modifed from 149, 123 and 102 MeV to 128, 111 and 107 
MeV respectively. It is seen that, at high densities there is an increase in the masses of 
both D~ and D" mesons in isospin asymmetric medium as compared to those in the isospin 
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symmetric nuclear matter for temperatures T=0, 50 and 100 MeV. However, at T = 150 
MeV, it is observed that for densities upto about 4.5po; the mass of D^ meson is higher in 
the isospin symmetric matter as compared to in the isospin asymmetric matter with rj = 
0.5. It is also seen that the modifications in the masses of D~ mesons is negligible as we 
change t] from to 0.3 upto a density of about 4po- For the Z)° meson, one sees that the 
isospin dependence is negligible upto 7^=0.3. This is because the drop in the mass of Z)° 
mesons due to isospin asymmetry given by Weinberg- Tomozawa term almost cancels with 
the increase due to the scalar and range terms as we go from rj = to rj = 0.3. At zero 



temperature 48[ as well as for temperatures T=50 and 100 MeV, there is seen to be an 
increase in the mass of the D mesons {D^ ,D^) as we go from isospin symmetric medium 
to the isospin asymmetric medium. This is because for T=0,50 and 100 MeV, the increase 
in mass of D given by the scalar exchange and the range terms dominate over the drop 
given by the Weinberg Tomozawa term as we go from symmetric nuclear medium (r^ = 0) to 
isospin asymmetric nuclear medium [r] =0.1,0.3 0.5). However, at T = 150 MeV, for rj =0.5, 
the drop given by Weinberg term dominates over the rise given by scalar and range terms 
for D^ and upto a density of about 4.5po for D~, and therefore mass of D meson decreases 
as we go from symmetric nuclear medium to isospin asymmetric nuclear medium in these 
density regimes. 

The medium modifications of the masses of D and D mesons in the present investigation 
are due to the interactions with the nucleons and scalar mesons a, ( and 6 in the hot 
nuclear medium. The values of the scalar fields in the medium are obtained by solving the 
the equations of motion for the scalar fields and the dilaton field x^ given by the coupled 
equations (9) to (12). The temperature and density dependence of the dilaton field x is seen 
to be negligible and thus the changes of the values of the scalar fields are observed to be 
marginal for the present investigation when the medium dependence of the x fi^ld is taken 
into account as compared to when its medium dependence is not taken into account (the 
so-called frozen glueball approximation where the value of x is taken to be its vacuum value). 
This leads to the modifications of the D and D meson masses as marginal as compared to 
the case when the medium dependence of the dilaton field is not taken into account. For 
example, for temperature T = 0, the drop in the mass of D+ mesons in isospin symmetric 



33 



nuclear medium is about 81 MeV and 364 MeV at ps = po and 4po respectively in the 
frozen glueball approximation, which may be compared to the values of mass drop as 78 
MeV and 347 MeV in the present investigation, when we take into account the effect of 
variation of the dilaton field with density. Hence one observes the difference between them 
to be marginal, of the order of about 5%. Similarly, the mass drop in the D~ meson for 
the isospin asymmetric matter at T=0 in the frozen glueball approximation is observed to 
be about 30 MeV and 184 MeV at pb = Po and 4po respectively. These are different from 
the values 27 MeV and 162 MeV of the present calculations with medium dependent dilaton 
field, by about 10%. 

In the present calculations, we have used the value of decay constant, fn = 135 MeV. 
In isospin symmetry nuclear medium, for temperature T =0, we observe that the values of 
mass drops for D~^ and D~ mesons are about 42 MeV and 4 MeV respectively at nuclear 



saturation density when we set /j:)=157 MeV [78|]. These values may be compared to the 
values of 78 MeV and 27 MeV when the fo is taken to be 135 MeV. Hence by modifying 
the value of fo by about 15% leads to modifications of the mass shifts of D^ and D^ by 
about 46% and 85% at density, po for symmetric nuclear matter at zero temperature. The 
mass drops for the value of the density as 4po for symmetric nuclear matter at T=0, are 
modified from 347 MeV and 162 MeV to the values 238 MeV and 93 MeV for the D+ and 
D^ respectively, when we change the value of fo from 135 MeV to 157 MeV. Hence, for 
Pb = 4po, the modifications for the D^{D^) masses are about 30% and 40% when we change 
the value of the D- meson decay constant. Hence, there seems to be appreciable dependence 
of the mass shifts of the D{D) mesons with the D-meson decay constant. 

We next examine how the masses of the D and D mesons change if we determine the 
parameters di and ^2 from DN scattering lengths calculated to be —0.43 fm and —0.41 fm 
in the 1 = and 1 = 1 channels respectively, in a coupled channel approach [56t]. In figures 
[9] and [lO], we show the variation of the energies of D and D mesons respectively, at zero 
momentum, with baryon density ps for different values of the isospin asymmetry parameter 
rj and for the values of the temperature as T = 0, 100, 150 MeV. The values of di and 
d2 parameters determined from these values of the DN scattering lengths turn out to be 
8.95/1711) and 0. 52/771 £> respectively, which can be expressed in terms of the mass of the kaon 
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FIG. 9: (Color online) The energies of D^ meson ((a),(c) and (e)) and of D^ meson ((b), (d) and 
(f)), at momentum k = 0, versus the baryon density (in units of nuclear saturation density), Pb/po, 
for different values of the isospin asymmetry parameter (ry = 0, 0.1, 0.3, 0.5) and for given values of 
temperature (T = 0, 100 MeV and 150 MeV). The values of parameters di and d2 are calculated 
from DN scattering lengths in I = and 1=1 channels. 
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FIG. 10: (Color online) The energies of D~ meson ((a),(c) and (e)) and of D^ meson ((b), (d) 
and (f)), at momentum /c = 0, versus the baryon density, expressed in units of nuclear saturation 
density, pb/ Pq, for different values of the isospin asymmetry parameter {rj = 0, 0.1, 0.3, 0.5) and for 
given values of temperature (T = 0, 100 MeV and 150 MeV). The values of parameters di and ^2 
are calculated from DN scattering lengths in 1 = and 1 = 1 channels. 
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as 2.385/mii- and O.lA/rriK- These values of parameters are smaller than the values of di and 
d2 as 2.56 /rriK and 0.73 /rriK respectively, when determined from the KN scattering lengths 
in J = and 1 = 1 channels. Since both of these terms are attractive, the masses of the 
D {D) mesons turn out to have a smaller drop when these parameters are fitted from DN 
scattering lengths as compared to when these are fitted from the KN scattering lengths. 

With the set of values of di and c?2 parameters as fitted from the DN scattering lengths, 
the drop in the masses oi D'^ meson in the isospin symmetric nuclear medium at nuclear 
matter saturation density po turns out to be 51, 42 and 40 MeV, at the values of temperature 
as, T = 0, 100 and 150 MeV respectively. These may be compared to the values of the mass 
drop of D^ meson of 78, 65 and 62 MeV for T=0, 100 and 150 MeV, respectively, when the 
parameters are fitted from the KN scattering lengths. For D^ mesons the values of mass 
drop are 50,42 and 40 MeV at T=0, 100 and 150 MeV respectively, when di and d2 are 
fitted from the DN scattering lengths. These may be compared to the values of mass drop of 
D^ meson of about 78, 65 and 62 MeV, when these are fitted to the KN scattering lengths. 
At higher densities, the di and d2 terms become more dominant and overcome the repulsive 
interaction of the first range term, leading to a drop of the D-meson masses due to the range 
term as well. At a density of ps = 4po, the drop in the D^ mass is observed to be about 286 
MeV, 256 MeV, 240 MeV for T=0, 100 MeV and 150 MeV, which may be compared to the 
values of 347 MeV, 313 MeV and 294 MeV for T=0,100, 150 MeV, when the parameters di 
and d2 are fitted from the KN scattering lengths. 

We plot the masses of the D mesons in figure [101 with the values of di and d2 are fitted 



from DN scattering lengths [56|. At the nuclear saturation density po, the value of mass of 
D^{D^) meson in isospin symmetric nuclear medium, is observed to increase by 2, 10 (10.5) 
and 12 (12.5) MeV at T = 0, 100 and 150 MeV, respectively. However, at higher densities, 
the di and d2 terms become more dominant thus leading to a drop of the D'{D^) masses 
in the nuclear matter. For ps = 4po, the mass of D^{D^) meson is seen to decrease by 87 
(86) MeV, 52 (51.5) MeV and 34 (33.6) MeV respectively. These may be compared to the 
results of the in-medium masses of D^{D^) meson, when di and d2 are fitted from the KN 
scattering lengths. In the latter case, as already mentioned for T=0, 100 MeV and 150 MeV 
and for ps = 4po, there is seen to be a drop of D^{D^) mass of 27.2 (27.2) 14.3 (14.2) and 
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10.8 (10.8) MeV at ps = po and 162 (162) 122.6 (122) and 101.7 (101.2) MeV at ps = 4po. 

As mentioned earlier, the in-medium mass of D^ meson decreases with increase in the 
isospin asymmetry of the nuclear medium. However, with values of di and ^2 parameters 
fitted from DN scattering lengths, the decrease in the mass with isospin asymmetry of the 
nuclear medium is observed to be smaller than that when these parameters are fitted from 
the KN scattering lengths. For example, at pb = 4po and T = 0, as we move from 77 = to 
1] = 0.3, the in-medium masses of D^ mesons decrease by 27 MeV with the values from the 
KN scattaring lengths and by 13 MeV with the set of values of di and ^2 parameters fitted 
from the DN scattering lengths. For the D^ meson, the di term gives a rise in the mass 
whereas the d2 term gives drop in the mass when one has a nonzero value of the asymmetry 
parameter as compared to the symmetric nuclear matter. Also, the first range term gives an 
increase due to isospin asymmetry, whereas both the Weinberg Tomozawa term and scalar 
meson exchange term give drop in the mass of the D^ meson. Due to the drop arising 
from the ^2 term, Weinberg Tomozawa term and the scalar term dominating over the the 
increase due to the di term, the mass of D^ meson decreases in the isospin asymmetric 
nuclear medium as compared to the mass in symmetric matter. When we use the values 
of the parameters di and ^2 as calculated from the DN scattering lengths, then there is a 
smaller contribution to the drop of the mass due to ^2 term, as compared to the drop due to 
isospin asymmetry arising to this term when the parameters di and ^2 are calculated from 
the KN scattering lengths. This is due to the smaller value for the ^2 in the former case. 
Therefore, with values of di and ^2 parameters calculated from the DN scattering lengths, 
the drop in the mass of D^ meson with isospin asymmetry of the medium is seen to be 
small. 

With the set of values of the parameters di and ^2 as fitted from the KN scattering 
lengths, the mass of the D~ meson is observed to increase with the isospin asymmetry of 
the medium. With the values of di and ^2 parameters as fitted from the DN scattering 
lengths, there is seen to be larger increase in the mass of the D~ meson due to the isospin 
asymmetry of the nuclear medium. For example, at Pb = 4po and T = 0, as we move from 
7] = to 1] = 0.3, the mass of the D^ meson is seen to increase by about 7 MeV when we 
use the values of di and ^2 as fitted from the KN scattering lengths and by 32 MeV with 
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the values of di and d2 fitted from the DN scattering lengths. For the D~ meson, there is 
an increase in the mass due to isospin asymmetry due to the Weinberg Tomozawa term, the 
first range term (term with coefficient (—7-) and the di term, whereas the scalar term and 
d2 term give drop in the mass of D~ meson as compared to the symmetric nuclear matter. 
The net effect is that the mass of the D~ meson increases with the isospin asymmetry of the 
nuclear medium. With the values of di and ^2 fitted from DN scattering lengths, because 
of the smaller value of (^2, the drop arising from the ^2 term due to isospin asymmetry in 
the medium is observed to be smaller than the case when the parameters di and ^2 are 
fitted from the KN scattering lengths (same as for D^ mesons). There is seen to be larger 
increase in the in-medium masses of D~ mesons as a function of the isospin asymmetry of 
the nuclear medium when the parameters are determined from the DN scattering lengths as 
compared to the KN scattering lengths. 

With the values of di and d2 calculated from KN scattering lengths, the mass of the D^ 
meson increases with the isospin asymmetry r] of the medium as shown in figure [HI However, 
if we use the values of di and d2-, fitted from the DN scattering lengths, then the mass of 
D^ mesons are seen to decrease with increase in the isospin asymmetry of the medium as 
shown in the figure [TOl For example, at pb = 4po and T = 0, as we move from 77 = to 
rj = 0.3, the D^ mass increases by about 10 MeV when the parameters di and ^2 are fitted 
from the KN scattering length, whereas the mass of D^ is seen to decrease by about 16 
MeV when di and d2 are calculated from the DN scattering lengths. The reason is, as a 
function of isospin asymmetry of the medium the di and d2 terms give rise to the masses of 
D^ mesons and the first range term gives a drop in the D^ mass. The values of di and d2 
are larger when fitted from the KN scattering lengths as compared to when calculated from 
the DN scattering lengths. Hence in the former case, the the increase in D^ from isospin 
asymmetry arising from the di and d2 terms dominates over the drop given by first range 
term. Therefore, with isospin asymmetry of the medium the mass of D^ increases in the 
former situation. However, for the values of di and d2 parameters fitted from DN scattering 
lengths, the increase due to di and ^2 terms is dominated by first range term. The Weinberg 
term gives a drop and the scalar term gives an increase in the mass of D'^ meson with the 
isospin asymmetry of the nuclear medium. The net effect on the mass of the D^ meson is a 
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drop with isospin asymmetry, rj, of the nuclear medium. 

The parameters di and d2 have the same effect on the masses of the D^ and D° mesons 
of giving an increase in their mass in the isospin asymmetric medium as compared to the 
masses in the symmetric nuclear matter whereas the first range term gives a drop in their 
masses. The Weinberg Tomozawa term and the scalar meson exchange term lead to an 
increase in the mass of D^ meson with isospin asymmetry. As can be seen from figures [7] 
and IHl the net effect on the D^ meson mass is an increase with isospin asymmetry, but the 
rise is less for the case when di and d2 are calculated from the DN scattering lengths due to 
the smaller values of di and d2 as compared to when these parameters are determined from 
the KN scattering lengths. For example, at ps = 4po and T = 0, as we move from 77 = to 
1] = 0.3, the in-medium mass of D^ mesons increases by about 45 MeV when the parameters 
are calculated from KN scattering lengths and by about 26 MeV when fitted from the DN 
scattering lengths. 

The mass modifications of D-mesons at finite density have been studied in the QCD sum 
rule and the mass shift at nuclear matter saturation density was found to be about —50 



MeV 37|. In the QMC model, the mass shift was around —60 MeV 39|]. In the present 
investigation, at finite densities the magnitude of the scalar fields first increases with increase 
in the temperature upto a particular value of temperature after which it starts decreasing. 
This behaviour is then reflected in the variation of the nucleon mass with temperature at 
finite densities. In QMC model, the behaviour of the scalar field is also seen to be the 

n 

same as in the present model [79[ . However, in the QMC model, the nucleon mass is seen to 
monotonically rise with temperature and there is no change in this trend observed even upto 
a temperature of about 250 MeV [79| unlike in the present chiral model or in the Walecka 
model 771]. This is because a field in the QMC model is not as strong as in the chiral model 
or the Walecka model. In QMC model, there are contributions to the masses of nucleons 
from the thermal excitations of the quarks inside the nucleon bag. This contribution of 



quarks dominates over the a field in QMC model 



79| . The small attractive mass shift 



for the D mesons, obtained within our calculations are in favor of charmed mesic nuclei 



as suggested in the QMC model [39|. This is, however, contrary to a repulsive potential 



obtained for the D mesons in the coupled channel approach 55| . In our investigation, if we 
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FIG. 11: (Color online) The optical potential of D^ meson ((a),(c) and (e)) and of D^ meson 
((b), (d) and (f)), are plotted as functions of momentum for ps = po, for different values of the 
isospin asymmetry parameter (77 = 0, 0.1, 0.3, 0.5) and for given values of temperature (T = 0, 100 
MeV and 150 MeV). The values of parameters di and ^2 are calculated from KN scattering lengths 
in I = and 1 = 1 channels. 
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FIG. 12: (Color online) The optical potential of D^ meson ((a),(c) and (e)) and of D^ meson 
((b), (d) and (f)), are plotted as functions of momentum for ps = po, for different values of the 
isospin asymmetry parameter (77 = 0, 0.1, 0.3, 0.5) and for given values of temperature (T = 0, 100 
MeV and 150 MeV). The values of parameters di and ^2 are calculated from KN scattering lengths 
in I = and 1 = 1 channels. 
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do not take into consideration the effect of the range terms on the in-medium properties of 
D mesons then at nuclear saturation density po and temperature T = 0, the mass of D^ and 
D^ mesons drop by 144 MeV and 96 MeV respectively in isospin symmetric nuclear medium 
[t] = 0). The large mass-shift shows the absence of repulsion due to total range term at 



nuclear saturation density. However, in coupled channel approach of Ref.J55|], a repulsive 
mass-shift of 11 MeV is given for D mesons. Figures [TT] and [T3] show the isospin dependence 
of the optical potentials for the D mesons as functions of the momentum, for densities 
Po and 4po respectively and for values of the temperature as T = 0, 100, 150 MeV. These 
optical potentials are plotted for the values of di and d2 calculated from the KN scattering 
lengths. Figures [T2] and IT^ illustrate the optical potentials for the D doublet. The isospin 
dependence of optical potentials is seen to be quite significant for high densities for the 
D-meson doublet {D^, D^) as compared to those for the D doublet. This is a reflection of 
the strong isospin dependence of the masses of the D-mesons as compared to the D as has 
been already illustrated in figures [7] and [HI For the D mesons, it is seen, from figure [HI that 
the masses of the the D~ meson and D^ meson for a fixed value of the isospin asymmetry 
parameter, rj are very similar, an observation which was seen earlier for the zero temperature 
case [4^. These are reflected in their optical potentials, plotted in figures [12] and [TH where 
one sees a maximum difference of about 5 MeV or so between D~ and D^ for ps = Po and 
about 10 - 15 MeV for pB = 4po. The present investigations of the optical potentials for 
the D and D mesons show a much stronger dependence of isospin asymmetry on the D 
meson doublet, as compared to that in the D meson doublet, as was already observed for 
the zero temperature case. However, when the parameters di and ^2 are fitted from the DN 
scattering lengths, then one sees a greater sensitivity to the isospin asymmetry in the masses 
of the D mesons as compared to the masses of the D mesons, as illustrated in figures [9] and 

m 

We shall now investigate how the behaviour of the dilaton field x in the hot asymmetric 
nuclear matter affects the in-medium masses of the charmonium states J /ip,ip (3686) and 
iIj{3770). In figures [151 [IS] and [TT] we show the shifts of the masses of charmonium states 
J /ip,ip (3686) and '?/'(3770) from their vacuum values, as functions of the baryon density for 
given values of temperature T and for different values of the isospin asjTiimetry parameter, 
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FIG. 13: (Color online) The optical potential of D^ meson ((a),(c) and (e)) and of D^ meson 
((b), (d) and (f)), are plotted as functions of momentum for ps = 4po, for different values of the 
isospin asymmetry parameter (77 = 0, 0.1, 0.3, 0.5) and for given values of temperature (T = 0, 100 
MeV and 150 MeV). The values of parameters di and ^2 are calculated from KN scattering lengths 
in I = and 1 = 1 channels. 
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FIG. 14: (Color online) The optical potential of D^ meson ((a),(c) and (e)) and of D^ meson 
((b), (d) and (f)), are plotted as functions of momentum for pB = 4po, for different values of the 
isospin asymmetry parameter (77 = 0, 0.1, 0.3, 0.5) and for given values of temperature (T = 0, 100 
MeV and 150 MeV). The values of parameters d\ and ^2 are calculated from KN scattering lengths 
in I = and 1 = 1 channels. 



45 



f]. We have shown the results for the values of the temperature, T = 0, 50, 100 and 150 
MeV. At the nuclear matter saturation density, ps = Po at temperature T = 0, the mass- 
shift for J/ip meson is —8.6 MeV in the isospin symmetric nuclear medium [r] = 0) and in 
the asymmetric nuclear medium, with isospin asymmetry parameter t] = 0.5, it is seen to 
be about —8.4 MeV. For ps = 4po and at zero temperature, the mass-shift for J/ip meson is 
observed to be about —32.2 MeV in the isospin symmetric nuclear medium [r] = 0) and in 
isospin asymmetric nuclear medium [r] = 0.5), it is seen to be modified to —29.2 MeV. The 
increase in the magnitude of the mass-shift, with density ps, is because of the large drop in 
the dilaton field x at higher densities. However, with increase in the isospin asymmetry of 
the medium the magnitude of the mass-shift decreases because the drop in the dilaton field 
X is less at a higher value of the isospin asymmetry parameter 77. For the nuclear matter 
saturation density ps = Po and at temperature T = 0, the mass-shift for ■?/'(3686) is observed 
to be about —117 and —114 MeV for rj = and 0.5 respectively, and for ■?/'(3770), the values 
of the mass-shift are seen to be about —155 MeV and —150 MeV respectively. At ps = 4po 
and zero temperature, the values of the mass-shift for ■i/'(3686) are modified to —436 MeV 
and —396 MeV for rj = and 0.5 respectively, and, for ■?/'(3770), the drop in the masses are 
about —577 MeV and —523 MeV respectively. As mentioned earlier, the drop in the dilaton 
field, Xj at finite temperature is less than at zero temperature and this behaviour is reflected 
in the smaller mass-shift of the charmonium states at finite temperatures as compared to 
zero temperature case. At nuclear matter saturation density ps = Po, temperature T = 100 
MeV, the values of the mass-shift for the J/ip meson are observed to be about —6.77 MeV 
and —6.81 MeV for isospin symmetric (77 = 0) and isospin asymmetric {1] = 0.5) nuclear 
medium respectively. At baryon density ps = 4po5 temperature T = 100 MeV, the mass- 
shift for J/ip is observed to be —28.4 MeV and —27.2 MeV for isospin symmetric {1] = 0) 
and isospin asymmetric {1] = 0.5) nuclear medium respectively. For the excited charmonium 
states ■?/'(3686) and ■?/'(3770), the mass-shifts at nuclear matter saturation density pb = Po and 
temperature T = 100 MeV, are observed to be —91.8 MeV and —121.4 MeV respectively, 
for isospin symmetric nuclear medium {1] = 0) and —92.4 MeV and —122 MeV for the 
isospin asymmetric nuclear medium with rj = 0.5. For a baryon density of Pb = 4poj and 
temperature T = 100 MeV, the mass-shifts for charmonium states ■?/'(3686) and (■?/'(3770)) are 
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seen to be —386 MeV and —510 MeV respectively, for isospin symmetric [r] = 0) and —369 
MeV and —488 MeV for isospin asymmetric nuclear medium with rj = 0.5. For temperature 
T = 150 MeV and at the nuclear matter saturation density Pb = Po, the mass-shifts for the 
charmonium states J /ip, 1^(3686) and iIj{3770) are seen to be —6.25, —85 and —112 MeV 
respectively in the isospin symmetric nuclear medium [i] = 0). These values are modified 
to —7.2, —98 and —129 MeV respectively, in the isospin asymmetric nuclear medium with 
T] = 0.5. At a baryon density pb = 4/9o, the values of the mass-shift for J/ip, ip{3686) 
and '?/'(3770) are observed to be —26.4, —358 and —473 MeV in isospin symmetric nuclear 
medium (77 = 0) and in isospin asymmetric nuclear medium with rj = 0.5, these values are 
modified to —27.6, —375 and —494 MeV respectively. Note that at high temperatures, e.g., 
at T = 150 MeV the mass-shift in isospin asymmetric nuclear medium {1] = 0.5) is more 
as compared to isospin symmetric nuclear medium [i] = 0). This is opposite to what is 
observed for the zero temperature case. The reason is that at high temperatures the dilaton 
field X has a larger drop in the isospin asymmetric nuclear medium [rj = 0.5) as compared 
to the isospin symmetric nuclear medium (77 = 0), due to the contributions from the 6 field 
for the nonzero rj, which is observed to decrease in its magnitude at high temperatures. 
The dependence of the wave function of the charmonium on the density and temperature 
of the medium can be introduced through the modification of the strength of the harmonic 
potential for the charmonium state [33| given as the parameter /3 in equation (30). In the 
medium, one expects the strength of the confining potential to be smaller than in the vacuum 
as due to medium modifications of the hadrons, more decay channels can become accessible, 
which are not available in the vacuum. We observe that when the parameter is decreased 
by about 5%, the mass drops of the charmonium states are increased by about 14% and a 
change of the parameter f3 by 2% leads to a mass drop of the charmonium states larger by 
about 5%. 

The values of the mass-shift for the charmonium states obtained within the present in- 
vestigation, at nuclear matter saturation density po and temperature T = 0, are in good 
agreement with the mass shift of J/ilj,tp{3686) and iIj{3770) as -8,-100 and —140 MeV 



respectively, at the nuclear matter saturation density computed in Ref. [3l| from the sec 



ond order stark effect, with the gluon condensate in the nuclear medium computed in the 
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FIG. 15: (Color online) The mass shift of J/ij^ plotted as a function of the baryon density in 
units of nuclear matter saturation density at given temperature, for different values of the isospin 
asymmetry parameter, rj. 

linear density approximation. On the other hand, in the present work, the temperature and 
density dependence of the gluon condensates are calculated from the medium modifications 
of the dilaton field within the chiral SU(3) model. In Ref. 3l|, the masses of the charmo- 
nium states were calculated for the symmetric nuclear matter at zero temperature, whereas 
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FIG. 16: (Color online) The mass shift of ^'(3686) plotted as a function of the baryon density in 
units of nuclear matter saturation density at given temperature, for different values of the isospin 
asymmetry parameter, r]. 

the present investigation studies the isospin asymmetry dependence of the masses of the 
charmonium states in the nuclear medium at finite temperatures, which will be relevant for 
the asymmetric heavy ion collision experiments planned at the future facility at GSI. The 
mass-shift for J/ip has also been studied with the QCD sum rules in 35| and the value at 
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FIG. 17: (Color online) The mass shift of ^'(3770) plotted as a function of the baryon density in 
units of nuclear matter saturation density at given temperature, for different values of the isospin 
asymmetry parameter, r]. 



nuclear saturation density was observed to be about —7 MeV. In [S^] the operator product 



expansion was carried out upto dimension six and mass shift for J/tp was calculated to be 
—4 MeV at nuclear matter saturation density po and at zero temperature. The effect of 
temperature on the J/-ip in deconfinement phase was studied in 75|, |80(] . In these investi- 
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gations, it was reported that J /ip mass remains essentially constant within a wide range of 
temperature and above a particular value of the temperature, T, there is seen to be a sharp 



change in the mass of J /if) in the deconfined phase. For example, in Ref. |8ll] the mass shift 
for J /ip was reported to be about 200 MeV at T = 1.05 Tc. In the present work, we have 
studied the effects of temperature, density and isospin asymmetry, on the mass modifica- 
tions of the charmonium states {J /ip-, ^(3686) and "0(3770)) in the confined hadronic phase, 
arising due to modifications of a scalar dilaton field which simulates the gluon condensates 
of QCD, within a chiral SU(3) model. The effect of temperature is found to be small for the 
charmonium states J/'?/'(3097), ^(3686) and ?/'(3770), whereas the masses of charmonium 
states are observed to vary considerably with density, in the present investigation. 

The medium modifications of the masses of D and D mesons as well as that of charmo- 
nium states could be an explanation for the observed J /ib suppression observed by NA50 
collaboration at 158 GeV/nucleon in the Pb-Pb collisions [20]. Due to the drop in the mass 
of the DD pair in the nuclear medium, it can become apossibility that the excited states of 



charmonium {ip , Xc2^ Xa, Xco) can decay to DD pairs [48| and hence the production of J /'^ 
from the decay of these excited states can be suppressed. Even at high values of densities at 
given temperatures, it can become a possibility that J /ip itself decays to DD pairs. Thus 
the medium modifications of the D mesons can modify the decay widths of the charmonium 
In figures [18] and [19], we show the density dependence of the masses of the 



states 



D^D^ as well as D^D^ pairs calculated in the present investigation, for temperatures, T = 
0, 100 and 150 MeV and for isospin asymmetry parameter, r^ = and 17=0. 5 respectively. 
We also show the in-medium masses of the charmonium states J/'?/', V^ (3686) and ip{?)77Qi) 
in these figures. We observe that in the isospin symmetric nuclear medium at zero tem- 
perature, the in-medium mass of charmonium ?/'(3770) is less than D^D^ and D^D^ pairs 
above baryon densities 0.6po and 0.8po respectively and therefore its decay does not seem 
possible at densities higher than these densities in the nuclear medium. However, as we 
move to the isospin asymmetric medium (77 = 0.5), the medium modifications for the masses 
of the DD pairs as well as of the charmonium states indicate that the decay of '?/'(3770) to 
D^D^ pairs can be possible above a density of about 2po, but the decay to D^D^ does not 
seem possible above a density of about nuclear matter saturation density. In the isospin 
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symmetric nuclear medium, the decay of the charmonium state tp{3Q86) to D^D^ and to 
D^D'^ seem possible above densities of about 3.4 po and 3.3 po respectively. In isospin asym- 
metric nuclear medium [i] = 0.5) the decay of the charmonium state iIj{3686) to D^D~ 
pairs seem as possibilities above a density of about 2po. The effects of the temperatures on 
the decay of the charmonium states to DD pairs have also been illustrated in the figures 
[18] and [191 The decay of '?/'(3770) to the DD pairs do not seem possible above a density of 
about Po even at T=100 and 150 MeV. However, for '?/'(3686), for T=100 MeV, the decay to 
the D^D^ and D^D^ seem possible above densities of about 4po and 3.8 po for symmetric 
nuclear matter For T=150 MeV, these values are modified to 4.5po and 4.3 po respectively, 
for 7^=0. As we move to the asymmetric nuclear matter, the densities above which the 
decay of ?/'(3686) decaying to D^D^ becomes possibile are 2.4 po and 2.2 po respectively. 
Similar to the zero temperature case, we do not see a possibility of '?/'(3686) to D^D^ for the 
asymmetric nuclear matter (7^=0. 5) at T=100 and 150 MeV. In the present investigation, 
the decay of J/^ to DD pairs does not seem as a possibilty even upto a density of about 
6po. We observe from figures 18 and 19 that the temperature dependence is minimal for 
the decay of the charmonium states to the DD pairs even though the density dependence 
is quite appreciable. 



The decay of the charmonium states have been studied in Ref. 29|, [33|. It is seen to 
depend sensitively on the relative momentum in the final state. These excited states might 
become narrow 33| though the D meson mass is decreased appreciably at high densities. 
It may even vanish at certain momentum corresponding to nodes in the wave function 33 |. 
Though the decay widths for these excited states can be modified by their wave functions, 
the partial decay width of Xc2i owing to absence of any nodes, can increase monotonically 
with the drop of the D^D^ pair mass in the medium. This can give rise to depletion in the 
J /^ yield in heavy-ion collisions. The dissociation of the quarkonium states (\E' ,Xc, J 1'^) 



into DD pairs has also been studied 82|, [83| by comparing their binding energies with the 



lattice results on the temperature dependence of the heavy-quark effective potential 84 1. 
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FIG. 18: (Color online) The masses of the DD pairs [D+D- in (a), (c), (e) and D^D^ in (b), (d), 
(f)] in MeV plotted as functions of ps/ Po for isospin for the symmetric nuclear matter (r/=0) and 
for temperatures, T = 0, 100, 150 MeV. 

VI. SUMMARY 

We have investigated in a chiral model the in- medium masses of the D, D mesons and the 
charmonium states {J /ip-, ^(3686) and '(/'(STTO)) in hot isospin asjTiimetric nuclear matter. 
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FIG. 19: (Color online) The masses of the DD pairs [D+D- in (a), (c), (e) and D^D° in (b), 
(d), (f)] in MeV plotted as functions of Pb/po for isospin asymmetry parameter value ij = 0.5 and 
temperatures, T = 0, 100, 150 MeV. 

The properties of the hght hadrons - as studied in SU{3) chiral model - modify the D{D) 
meson properties in the dense and hot hadronic matter. The SU{?)) model, with parameters 
fixed from the properties of the hadron masses in vacuum and low-energy KN scattering 
data, is extended to SU(4) to derive the interactions of D{D) mesons with the light hadron 
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sector. The mass modifications of D^ and D^ mesons is strongly dependent on isospin- 
asymmetry of medium when we determine the parameters di and d2 consistent with the KN 
scattering lengths. However, the sensitivity to the isospin asymmetry is seen to be more 
for the D doublet, when we fit the parameters to the DN scattering lengths as calculated 
in the coupled channel approach in Ref. 56|. At finite densities, the masses of D{D) 
mesons are observed to increase with temperature [47| upto a temperature above which it 



is observed to decrease. The mass modification for t 
earlier finite density calculations of QCD sum rules 



le D mesons are seen to be similar to 
38l . l74l | as well as to the quark-meson 



coupling model |39| 



This is in contrast to the small mass modifications in the coupled 



channel approach |5lL l55| . Also we obtain small attractive mass shifts for D mesons similar 



to the results obtained from the QMC model, which might lead to formation of charmed 
mesic nuclei. These results for the D mesons are contrary to the results from the coupled 
channel approach 55|] , where the D mesons experience a repulsive interaction in the nuclear 
medium. In our calculations the presence of the repulsive first range term (with coefficient 
— J- in Eq. ( TT8|) ) is compensated by the attractive di and ^2 terms in Eq.( TT8i) . Among the 
attractive range terms {di and d2 terms), di term is found to be dominating over d2 term. 

We have investigated in the present work, the effects of density, temperature and isospin 
asymmetry of the nuclear medium on the masses of the charmonium states J/^, '?/'(3686) 
and (^^(3770), arising due to modification of the scalar dilaton field, Xi which simulates 
the gluon condensates of QCD, within the chiral SU(3) model. The change in the mass of 
J /ip with the density is observed to be small at nuclear matter saturation density and is 
in agreement with the QCD sum rule calculations. There is seen to be appreciable drop 
in the in-medium masses of excited charmonium states ?/'(3686) and '0(3770) with density. 
The mass drop of the excited charmonium states '?/'(3686) and '?/'(3770) are large enough 
to be seen in the dilepton spectra emitted from their decays in experiments involving p-A 
annihilation in the future facility at GSI, provided these states decay inside the nucleus. 
The life time of the J /ip has been shown to be almost constant in the nuclear medium, 
whereas for these excited charmonium states, the life times are shown to reduce to less than 



5 fm/c, due to appreciable increase in their decay widths 33(|. Hence a significant fraction of 
the produced excited charmonium states in these experiments are expected to decay inside 
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the nucleus [85|. The in-medium properties of the excited charmonium states '?/'(3686) and 
iIj{3770) can be studied in the dilepton spectra in p-A experiments in the future facihty of 
the FAIR, GSI [16]. The mass shift of the charmonium states in the hot nuclear medium 
seem to be appreciable at high densities as compared to the temperature effects on these 
masses, and these should show in observables like the production of these charmonium states 
in the compressed baryonic matter experiment at the future facility at GSI, where baryonic 
matter at high densities and moderate temperatures will be produced. 

The medium modifications of the D meson masses can lead to a suppression in the J/\E' 
yield in heavy-ion collisions, since the excited states of the J/"^ can decay to DD pairs in the 
dense hadronic medium. The medium modifications of the masses of the charmonium states 
as well as the D and D mesons have been considered in the present investigation. The isospin 
asymmetry lowers the density at which decay to D^D^ pairs occur. Due to increase in the 
mass of D^D^ in the isospin-asymmetric medium, isospin-asymmetry is seen to disfavor 
the decay of the charmonium states to the D^D^ pairs. At zero or finite temperatures, 
there does not seem to be a possibility of decay of J/\E' to D^D^ or D^D^ pairs. The 
isospin dependence of D~^ and D^ masses is seen to be a dominant medium effect at high 
densities, which might show in their production {D^/D'^), whereas, for the D~ and D^, 
one sees that, even though these have a strong density dependence, their in-medium masses 
remain similar at a given value for the isospin-asymmetry parameter rj. This is the case 
when we fit the parameters di and d2 from the KN scattering lengths. When we determine 
these parameters from the DN scattering lengths as calculated in Ref. j56|], the masses of 
the D doublet are seen to be more sensitive to isospin asymmetry in the medium. The 
strong density dependence as well as the isospin dependence of the D{D) meson optical 
potentials in asymmetric nuclear matter can be tested in the asymmetric heavy-ion collision 
experiments at future GSI facility [16| in observables like the D^ / D^ as well as D^ / D^ 
ratios. In the present work, we have investigated the in-medium masses of the charmonium 
states due to their interaction with the scalar dilaton field (simulating the gluon condensates 
of QCD) and D{D) mesons due to the interaction with nucleons as well as scalar mesons. 
The parameters of the asymmetric nuclear matter are fitted to the nuclear matter properties, 
vacuum baryon masses, the hyperon potentials as well as the KN scattering lengths in the 
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chiral SU(3) model. The study of the in-medium modifications of D mesons in hadronic 
matter including hyperons along with nucleons at zero and finite temperatures, as well as 
the study of the medium modifications of the strange open charm mesons can be possible 
extensions of the present investigation. 
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